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The influence of electron-phonon interaction on magnetotransport in two-dimensional electron 
systems under microwave irradiation is studied theoretically. Apart from the phonon-induced re- 
sistance oscillations which exist in the absence of microwaves, the magnetoresistance of irradiated 
samples contains oscillating contributions due to electron scattering on both impurities and acoustic 
phonons. The contributions due to electron-phonon scattering are described as a result of the in- 
terference of phonon-induced and microwave-induced resistance oscillations. In addition, microwave 
heating of electrons leads to a special kind of phonon-induced oscillations. The relative strength of 
different contributions and their dependence on parameters are discussed. The interplay of numer- 
ous oscillating contributions suggests a peculiar magnetoresistance picture in high-mobility layers 
at the temperatures when electron-phonon scattering becomes important. 
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I. INTRODUCTION 

Studies of high-mobility two-dimensional (2D) electron 
systems in weak perpendicular magnetic fields B reveal 
several kinds of magnetoresistance oscillations which are 
caused by the scattering-assisted electron transitions be- 
tween different Landau levels and survive increase in 
temperature T, in contrast to Shubnikov-de Haas os- 
cillations. Most of attention is paid to the microwave- 
induced resistance oscillations (MIRO) which appear un- 
der steady-state microwave (MW) irradiation of 2D elec- 
tron gas.^ The MIRO periodicity is determined by the 
ratio of the radiation frequency uj to the cyclotron fre- 
quency Wc, while the amplitude is governed by the ra- 
diation power. As the power increases, MIRO minima 
are spread into intervals of zero dissipative resistance 
known as "zero-resistance states" .^''^'^ The basic theo- 
retical consideration^"^ of MIRO involves the picture of 
transitions between broadened Landau levels under elas- 
tic scattering of electrons by impurities or other inhomo- 
geneities in the presence of MW excitation, when elec- 
trons absorb or emit radiation quanta. As a consequence 
of these transitions, both the electron scattering rate 
and the electron distribution function acquire oscillating 
components which lead to magnetoresistance oscillations. 
The corresponding contributions are often described in 
terms of " displacement" ^^^'^ and "inelastic"®'^ micro- 
scopic mechanisms of MIRO, respectively. The above 
terminology emphasizes that the scattering-assisted MW 
absorption gives rise to spatial displacement of electrons 
along the applied dc field, while the modification of the 
distribution function is controlled by inelastic relaxation 
of electrons. Consideration of these mechanisms satis- 
factory explains both the periodicity and the phase of 
MIRO observed in experiments, as well as their power 
and temperature dependence. ^"'^^'^^ 

Inelastic scattering of electrons by acoustic phonons 
makes possible the phonon-assisted transitions of elec- 
trons between Landau levels. The probability of such 



transitions is maximal when the phonon momentum is 
close to the Fermi circle diameter, 2pF- This leads to 
a special kind of magnetophonon oscillations also known 
as phonon-induced resistance oscillations (PIRO),^^~^^ 
whose periodicity is determined by commensurability 
of the characteristic phonon energy with cyclotron en- 
ergy. In GaAs quantum wells with very high mobility 
(^ 10^ cm^/V s), these oscillations are visible already at 
T = 2 K,^^ while in the samples with moderate mobilities 
(- 10'5 cmVV s) they are weU seen at T = 10 - 20 K.^^ 
A theoretical description of PIRO assuming interaction 
of 2D electrons with anisotropic bulk phonon modes and 
showing a good agreement with experimental data has 
been presented recently. -"^^ 

Both MIRO and PIRO are observed in the same in- 
terval of magnetic fields, i? < 1 T, and have comparable 
periods. Therefore, it is interesting to study the inter- 
play of these remarkable oscillating phenomena. In other 
words, one may pose a question: what happens with 
magnetoresistance of the sample which shows PIRO if 
this sample is irradiated by microwaves? More gener- 
ally, there exists an important problem of the infiuence 
of electron-phonon interaction on MW-induced magne- 
toresistance. Some aspects of this problem have been 
considered previously. First, it is well established that 
increasing MW power leads to heating of the electron 
system, especially in the samples with lower mobility. 
The effective electron temperature is determined by 
the power transmitted from electron system to the lattice 
via electron-phonon interaction.-'^^ The heating causes a 
suppression of MIRO amplitudes"'2°'2i foj. ^-^^ 

mam rea- 
sons: the enhancement of Landau level broadening and 
a decrease of the inelastic scattering time. In this sense, 
a consideration of electron-phonon interaction becomes 
necessary for description of MW photoresistance at ele- 
vated MW power. Further, as shown theoretically,^^'^'^ 
there is a direct infiuence of electron-phonon scattering 
on the photoresistance: this scattering in the presence of 
MW excitation can give rise to pronounced magnetore- 
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sistancc oscillations whose behavior is different from that 
of the case of impurity-assisted transport. A considera- 
tion of clectron-phonon interaction as a possible reason 
for the absolute negative resistivity and zero-resistance 
states has been also presented. 

The aim of this paper is to give a consistent theoreti- 
cal description of magnetotransport under MW irradia- 
tion in the presence of both elastic and phonon-assisted 
electron scattering. It is (lemonstrated that the magne- 
toresistance has contributions of MIRO, PIRO, and the 
terms corresponding to interference of these types of os- 
cillations. Apart from this, heating of electron system 
leads to an additional PIRO contribution which is shifted 
by phase with respect to equilibrium PIRO. Superposi- 
tion of all these contributions gives rise to a complicated 
oscillating resistivity picture which is essentially differ- 
ent from either equilibrium magnetoresistance or MW- 
modified magnetoresistance under elastic scattering only. 
The analytical expression for magnetoresistance is ob- 
tained by assuming weak excitation, when the response 
is linear in MW power. 

The paper is organized as follows. Section II describes 
the quantum Boltzmann equation for the system of 2D 
electrons interacting with impurities and phonons in the 
presence of microwaves. Section III is devoted to cal- 
culation of electron distribution function and resistivity. 
Section IV contains numerical results based on the ex- 
pressions derived in Sec. Ill, description of different con- 
tributions to oscillating magnetoresistance, a discussion 
of the approximations used, and a brief conclusion. 



II. GENERAL FORMALISM 

In the region of weak magnetic fields, when the cy- 
clotron energy is small in comparison with the Fermi en- 
ergy, it is convenient to calculate the resistivity by us- 
ing the quantum Boltzmann equation for the generalized 
Wigner distribution function fept, which depends on en- 
ergy e, quasiclassical 2D momentum p and time t. The 
coordinate dependence of this function is omitted, since 
only spatially-homogeneous 2D systems are under con- 
sideration. To take into account the influence of MW 
radiation and dc excitation on 2D electrons, one may use 
a transition to the moving coordinate frame (see Ref. 9 
and references therein), when the MW and dc field po- 
tentials are transferred to the collision integrals and do 
not appear in the left-hand side of the kinetic equation. 
The scattering of electrons by impurities and phonons 
is treated within the self-consistent Born approximation 
(SCBA), which assumes that both the mean free path 
of electrons and the magnetic length are larger than the 
correlation lengths of the corresponding scattering po- 
tentials. Then, the quantum Boltzmann equation takes 
the form 
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where J*™, Jf^, and J^'^ are the collision integrals for 
electron-impurity, electron-phonon, and electron-electron 
scattering, respectively. They are written as 



dre'- 



dt' 
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where k is the index of corresponding interaction (im, 
ph, or ee). Throughout the paper, the system of units 
with fi, = 1 is used. The distribution function in the 
double-time representation, /ti.taCp); related to the 
Wigner distribution function f^-pt by the Wigner trans- 
formation, which is defined for arbitrary function A 
as Aept = / (ire*^'^A(_,_T-/2,t-r/2(p)- The Keldysh self- 
energies are described by the following expressions: 

S^:';^+(P) = - / (^"^(9)e^''-^*^*^G,-r*.(p-q), (3) 
X j ^2),-+(n)e-^"(*-*^)G,-r,.(p-q),(4) 



and 



s^:t+(p)=-2/(^t7x.rt.(p+q) 

X /(^G,-+,(p'-q)G+:,,(p')- (5) 

The retarded and advanced self-energies, S^**)^ and 
E^**)"*, can be expressed as 



(6) 



where 



S^:'^r^(p) = J^w{q)e^^-^^^^^Gt^%{p-^), (7) 
X J ^iP±±(fi)e-^(*-*^)G±±,(p-q), (8) 



and 



(9) 
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The quantities entering Eqs. (3)- (9) arc described be- 
low. First, G^^^(p) is the Keldysh Green's function for 
electrons, which is related to fti,t2{p) as 

Gt,%{p) = J dt' (p)G^,*,(p) - Gt1,t-(p)/t',t.(p)] , 

(10) 

where G and are the retarded and advanced Green's 
functions. The Green's function G"' is expressed by 
the linear combination, G"' = G + G^ — G^. Next, 
ZJ^Q (f2) is the Keldysh Green's function for phonons: 



+{N^^^ + l)6{n + ujxQt)], 



(11) 



where N^^^^ = [e^^^/-^ — is the Planck's distribution 
function and io\q is the phonon frequency. To find G , 
G++, D , and D^^ , one uses 



G~=G-++G^, G++=G-+-G^, 
D—=D-+ + D^, D++ = D-+-D^. 



(12) 



The retarded and advanced Green's function of phonons 
are expressed as 



D 



■ waq + iO fl + LUxQ + iO 



(13) 



and D^f^{Q) = [D^^^iQ)]* . The retarded Green's func- 
tion of electrons depends on the magnetic field B and 
satisfies the equation 



s — 
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2m Smc^ dp"^ 
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G%, = l + \jdTe'^^jdif (14) 
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where m is the effective mass of electrons. The equation 
for the advanced Green's function G^ differs from Eq. 
(14) only by a replacement of the superscripts R with A. 
The integral on the right-hand side of Eq. (14) describes 
the influence of impurity scattering on the energy spec- 
trum of 2D electrons in the magnetic field (the influence 
of other kinds of interaction on the electron spectrum is 
neglected here). 

The electron-impurity interaction in Eqs. (3) and (7) 
is described by the Fourier transform of the random po- 
tential correlator, w{q). The interaction with phonons 
[Eqs. (4) and (8)] is considered under approximations 
of equilibrium phonon distribution and bulk phonon 
modes. The influence of electron-phonon interaction on 
the phonon spectrum is neglected. The phonons are char- 
acterized by the mode index A and three-dimensional 
phonon momentum Q — (q, (7z). The squared matrix 
element of electron-phonon interaction potential is repre- 
sented as Maq = Gaq/(9z), where I{q,) = \ (0|e*9»^|0) 1^ 
is determined by the confinement potential which defines 
the ground state of 2D electrons, |0), and the function 



Gaq is determined by both deformation-potential and 
piezoelectric mechanisms of interaction: 



GxQ 
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Here V is the deformation potential constant, hu is the 
piezoelectric coupling constant, and Pm is the material 
density. The sums are taken over Cartesian coordinate 
indices, e^Qi are the components of the unit vector of 
the mode polarization, and the; coefficient is equal 
to unity if all the indices i,j,k are different and equal 
to zero otherwise. The polarization vectors and the cor- 
responding phonon mode frequencies are found from the 
eigenstate problem 



(16) 



where Kij (Q) is the dynamical matrix and 5ij is the Kro- 
necker symbol. For cubic crystals and in the elastic ap- 
proximation, the dynamical matrix is expressed through 
three elastic constants which are written in the conven- 
tional notations: 

Kij{Ci) = [(Cll - Caa)QI + CAAQ^]5i3 

+{ci2 + Cii)QiQj{l-5ij). (17) 

The 2D vector Rtjts standing in the exponential fac- 
tors in Eqs. (3), (4), (7), and (8) includes both the 
term with dc field E = {E^^Ey) and the MW- induced 
term which depends on the frequency w and amplitude 
E^ of the MW field strength. It is convenient to use 
= Rx ±iRy given by 



eE+ih-t2) eE^ 



mwc 



-iuiti 



— e 



and R 



tlt2 



where E± = E, 



± iEy. The form 



of the complex coefficients s± depends on polarization of 
the radiation. In the case of linear polarization (used for 
the calculations below). 



1 



(19) 



while in the case of circular (— ) polarization s+ = and 
S- = 1/ {u! — ujc + ii-Op)- Here ujp is the radiative de- 
cay rate which determines the cyclotron line broadening 
because of electrodynamic screening effect. ^^'^^'^^ Under 
usual experimental conditions, when the electromagnetic 
radiation is normally incident from the vacuum on the 
sample with dielectric permittivity e, and 2D electron 
layer is close to the surface of the sample, one can find 
Up = 2Tre^ns/mc\/e*, where ^/e* = (1 -|- •\/e)/2 and Ug is 
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the sheet electron density. It is assumed that iVp is much 
larger than the transport scattering rate of electrons, this 
condition is amply satisfied for high-mobility electrons. 

In contrast to J*'" and Jp'^ described above, the 
electron-electron collision integral J^^ is not affected by 
the external fields. In Eqs. (5) and (9) the "exchange" 
Coulomb terms in the self-energies are neglected, since 
the main contribution to the electron-electron collision 
integral at low temperatures comes from small-angle scat- 
tering. The dynamical screening effects are also ne- 
glected. The statically-screened 2D interaction potential 
is written as 



2776^ 



e(9 + qo) ' 



(20) 



where qo = 2e^m/e is the inverse screening length. 

Having solved the kinetic equation (1), one can find 
the dissipativc current density j , averaged over the period 
tu = 2Tr/uj of the microwaves; 



. _ 2e f de 



dp 

(2^ 



1 



dtG-+. 



(21) 



The current can be also expressed through the distribu- 
tion function /^pj with the aid of Eq. (10). The factor 
of 2 in Eq. (21) accounts for the spin degeneracy of elec- 
trons (the Zeeman splitting is neglected). 



III. CALCULATION OF RESISTIVITY 

It is assumed in the following that the temperature T 
is small in comparison with the Fermi energy ep- There- 
fore, all the quantities whose dependence on the abso- 
lute value of electron momentum is slow can be taken 
at p = |p| = Pf, where pp is the Fermi momentum re- 
lated to the sheet electron density rig as pp = \/2nns- 
Accordingly, instead of the function f^pt, one may use 
the distribution function f^^t which is equal to f^pt at 
p = Pf and depends on the angle (f of the quasiclassical 
momentum in the 2D plane. The absolute value of the 
momentum q = p — p' transferred in the elastic scat- 
tering by impurities is then given hy q = 2pFsiii[9/2), 
where 6 = (p — <p' is the scattering angle. In a similar 
way, since the characteristic phonon energy ujxq is small 
compared to the Fermi energy, the electron-phonon scat- 
tering is treated in the quasielastic approximation, with 

Q = \J qi+ 'ip% sin^ (0/2). The angle of the vector q is 
given by (fq = 7r/2 -t- (j), where (j) = {ip + ip')/2. 

Another approximation used below is the neglect of 
temporal harmonics of the distribution function fe^t- 
Strictly speaking, excitation of such harmonics by mi- 
crowaves can contribute to the dc current through 
the processes of second order in interaction, which 
leads to the so-called " photovoltaic" mechanism of MW 
photoresistance.^ According to theoretical estimates, this 
contribution is not large, and since no experimental evi- 
dence for the importance of the photovoltaic mechanism 



has been found, its neglect is justified. In the approxima- 
tions described above, Eq. (1) is reduced to a stationary 
kinetic equation for time-independent distribution func- 
tion /e^: 



, , ^/ey Tim 
i^c— — 



dip 



jph 

•-'eip 



j: 



sip- 



(22) 



The collision integrals standing here are calculated by 
using the equations of Sec. II (for transformation of 

electron-impurity collision integral, see also Refs. 9 and 
27). This leads to the following expressions: 



Jim 



2- d^' 



27r 

-^s-\-noj-\--y[fs-\-noj-\--y (p' /ey?] ; 



(23) 



X,n 

X {[(^-.Q 

~-^u:xQfsip]De+ujxQ+nuj+t} ■ (24) 



+ feip)fe-ujxQ+nu)+j ip' 



The rate i'{6) characterizing the impurity scattering is 
defined as v{9) = miuppi? sin(6'/2)], and the quantities 
LO\Q and MxQ are expressed through the variables qz, 0, 
and (j) as described above. Next, Jn{P) is the Bessel func- 
tion, and Ds = m~^ J j^{Gfp - Gfp)/i is the density 
of electron states expressed in units m/7r (so = 1 at 
B = Q). In the presence of magnetic field (wc ^ £f is 
assumed), the SCBA leads to the form 



Ds = 1 + 2 ak cos 



27rfc£ 



(25) 



k=l 



dk = ("1)*^^ ^ e^p{—Trk/L0cTq)Lk_i{2'Kk/uJcTc 



97' 



where is the quantum lifetime of electrons and L are 

the Laguerre polynomials. In the following, the case of 
overlapping Landau levels corresponding to weak mag- 
netic fields is considered. The density of states is then 
given by the expression = 1 — 2dcos{2TTe/u)c), where 
d = exp(— Tr/wcTg) is the Dingle factor, and Tq is the 
quantum lifetime of electrons due to impiirity scatter- 
ing: T^^ =■ v(9) (here and below, the line over the func- 
tion denotes the angular averaging (27r)^^ dO...). The 
angular-dependent variables /? and 7 standing in Eqs. 
(23) and (24) are proportional to the MW and dc fields, 
respectively: 



/3 = 



eE^VF 
\/2w 



s_ (e*v - e*^ ) -h s+ (e-*^ - e"'^ ) (26) 



and 



7 = 



evp 
2iijJc 



E_ (e*'^ - e'-^ ) - E+ie-'"^ - e-'"^ )], (27) 
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where vp = PF/m is the Fermi velocity. 

It is convenient to expand the distribution function in 
the angular harmonics: /g^ = X]fe /efc^*'^'^. The density 
of dissipative electric current is determined by the fc = 1 
harmonic. Applying the usual notation j± = ± ijy, 
one has 



IT ' 



(28) 



In the regime of classically strong magnetic fields, the 
anisotropic (A; ^ 0) part of the distribution function is ex- 
pressed through the isotropic (angular-independent) part 
fe = ,/e.fe=o directly from Eq. (22), by neglecting the an- 
gular dependence of all distribution functions in the col- 
lision integrals (23) and (24). The collision integral J®^ 
does not lead to angular relaxation of the anisotropic 
distribution and, therefore, does not contribute to such a 
solution. A subsequent linearization of the collision inte- 
grals (23) and (24) with respect to the dc field E allows 
one to find a linear response to this field in Eq. (28). 
The response to E- defines the symmetric part of the 
diagonal conductivity, j_ = aaE-, where 



(29) 



with 



^en — I^eDe+nuj-^ + [/e — fe+nu)] D^+nu] (30) 



and 



dfe 



+ [(^^c..Q + l)/e-(A^c..Q+/e)/e 



■ de 



D. 



e-uixQ+nuj 



.(31) 



The resistivity p is expressed as p ~ (t^ /cr^ , where a± = 
e^ns/mujc is the classical Hall conductivity. 

In the absence of MW radiation, when (3 = 0, one 
should substitute in Eq. (29) [J„(/3)]^ = 1 for n = 0, 
[J„(/3)]^ = for n 7^ 0, and take into account that is 
the equilibrium Fermi distribution function. Then, \is- 
ing Rl^ = Dlidfjde), and R^J = -(iV„,^ + 1)A(1 - 
fe-uxci)DeDs-i^^^/T, one obtains the expression for the 
linear resistivity of 2D systems. The phonon-assisted 
contribution, which is proportional to the factor iJ^Q , in 
this case is given by Eqs. (8) and (9) of Ref. 18. 

The MW excitation leads to a significant modification 
of the resistivity. First, one should consider in Eq. (29) 
the terms with nonzero n, corresponding to absorption 
of n quanta of radiation. Next, it is necessary to account 
for the changes of electron distribution function /g owing 



to this absorption. This modified distribution function is 
found from the isotropic kinetic equation 



jira ^ jph ^ jee ^ q_ 



(32) 



The collision integrals J*"* and Jf'* are given by Eqs. 
(23) and (24) , where 7 = and only the isotropic part of 
the distribution function is retained (since the anisotropic 
part is small, /e,^' ~ fs,p — fe). The isotropic electron- 
electron collision integral is given by®'^^ 



de' / dn A,,,in) [(1 - /,)(! - /eO,/s+n/.'-n 



-/e/.'(l - /e+n)(l - /e'-n)] D,,D,+aD, 



(33) 



and the function Aee'i^) is estimated as A ~ 
(2776^)^^ ln(i5'o/(Zi)5 where qi is a characteristic (small) 
momentum separating the ballistic and diffusive regimes 
of electron-electron scattering; see Ref. 8 for a more de- 
tailed consideration. This energy-independent form of A 
is valid in the experimentally relevant region of temper- 
atures, where the transferred momentum q ~ T/vp lies 
in the interval go ^ 9 ^ 91 ■ 

Below, only the terms linear in MW power are taken 
into account. Following the basic idea of Ref. 8, the 
distribution function is presented in the form 



/e - ft + 5fe 



(34) 



where ff slowly varies with energy on the scale of Wc, 
and (5/e oscillates with the period lOc- The function 
differs from the equilibrium Fermi distribution because of 
the influence of microwaves and satisfies Eq. (32) where 
the oscillating contribution to the density of states is ne- 
glected, Dg = 1. 

One should notice that the electron-electron scattering 
is stronger than electron-phonon scattering and controls 
the electron distribution if the MW intensity is not high. 
Thus, the function f^ can be reasonably approximated 
by the heated Fermi distribution characterized by the 
electron temperature . This temperature is determined 
from the balance equation 



(35) 



where Wq is the MW power absorbed by the electron sys- 
tem and Wr is the power transmitted out of this system 
through the energy relaxation owing to electron-phonon 
scattering. The balance equation can be obtained by 
integrating the kinetic equation (32), multiplied by the 
energy and density of states, over energy e. This leads 
to the following expressions: 



.tr 
im 



(36) 



where z/*^ = iy{0){l — cos0) is the transport rate of elec- 
trons due to electron-impurity scattering and 



(37) 
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is the dimensionless function proportional to MW power. 

The rate 



2 V^^'a;;,Q 



.)} (38) 



characterizes energy absorption due to electron-phonon 
scattering. The integral operators <S„ are defined as 



xrnMAQ(l -cos6')"A (39) 

where A is an arbitrary function which depends on the 
variables of integration. The function = [e^/"^" — 
difii'ers from the Planck's distribution function Nq 
by substitution of electron temperature in place of the 
lattice temperature T. The introduction of the operators 
<S„ considerably simplifies notations in the following. The 
power absorbed by the lattice is 



Wr 



m 



(40) 



Assuming Tg/T—l <C 1, one gets Wr oc T^—T. According 
to the balance equation (35), this means that the heating 
is directly proportional to the MW power. Since only the 
contributions linear in power are considered, the relation 
Te/T — 1 ^ 1 must be satisfied. In this approximation, 
the balance equation is rewritten as 



T 

Y 



\2TJ 



2 u^' 



^p\(w) 



(41) 



where the energy relaxation rate u^f^ is introduced by the 
expression 



(42) 



and F{x) = [a;/sinh(a;)]2. The function F{loxq/^T) 
can be also expressed through the Planck's distribution, 

F{ojxci/2T) = {u^xQ/T)^N^,ciN^,c + !)• In the high- 
temperature limit, when 2Te ^ WjWaq, the dependence 
Wr (X Te — T takes place even at high MW power. In 
this limit, //"^(cj) becomes frequency-independent and co- 
incides with the transport rate due to electron-phonon 
scattering. This rate is defined as 



(43) 



Once the distribution function is known, one can 
find the rapidly oscillating correction Sf^ from the os- 
cillating (proportional to the Dingle factors) part of Eq. 



(32). This leads to the following expression 



Sfe 



x{N^,^ - AT^^^+J - 2g^,^{N^,^ - N^^^)} , (44) 
where 



e-\-uj 



.m, (45) 



and 6Dg = — 1 = —2dcos{2'iTe/ujc)- The inelastic 
scattering time Ti„, which describes relaxation of the 
isotropic part of the distribution function due to both 
electron-electron and electron-phonon scattering, is de- 
fined as I /Tin = l/'^ee + ^/Tph- Thc clcctron-electron 
scattering contribution has been widely discussed in the 
literature^'^'"'^^. The result for the electron-phonon scat- 
tering contribution is 



-L^So{2N^,^ + l + f^^^^^ 

'ph 



/I-^aJ- (46) 



In derivation of this expression, the integrals 
So{{2Ncj^^ + 1) exp(±27rcjAQ/t^c)} containing rapidly 
oscillating factors are neglected, which allows one to 
consider only the outscattering contribution in the 
collision integral. At high temperatures, 2T > uj\q, 
energy dependence of Tp/j becomes inessential. For 
typical parameters of GaAs quantum wells, 1/Tph is 
small in comparison to l/ree, so Tin — Tee- According to 
Eq. (44), the function 5f^ is composed from three terms: 
the first one is caused by electron-impurity scattering, 
while the two other terms are due to electron-phonon 
scattering. The second term does not oscillate with 
MW frequency and exists owing to electron heating. 
The third term is directly proportional to MW power 
and, in this sense, is analogous to the first term. Since 
only the effects linear in MW power are considered, the 
saturation^'^ of 5f^ is neglected. Strictly speaking, this 
approximation also requires a neglect of heating in the 
third term on the right-hand side of Eq. (44), because 
accounting for ^ T would produce the contributions 
of higher orders in MW power in this term. 

Having found the distribution function, one can cal- 
culate the resistivity from Eqs. (29)-(31). It is written 
as 



/9 = Pe + Sp, 



(47) 



where pe and Sp are caused, respectively, by the compo- 
nents of the distribution function and Sf^. Consider 
first thc modifications of resistivity associated with the 
fimction f^. The heated Fermi distribution is substi- 
tuted into Eqs. (30) and (31), and only the terms linear 
in MW power are retained (n = 0, ±1). The integral over 
e in Eq. (29) is calculated under the approximation 



2w^T » Wc, 



(48) 
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when the Shubnikov-de Haas oscihations are thermally 
averaged out. The result of this integration can be di- 
vided into 5 parts: 

Pe = pL + pIh + pO^ + Pph + Pfh ■ (49) 

Two of these terms come from impurity-assisted (index 
im) scattering, the first one is the resistivity without ir- 
radiation, and another one is a well-known^'^'^^^^* MW- 
induced correction which describes MIRO due to dis- 
placement mechanism of photoresistance: 



Pim 



„dis 
Pim 



2 TTUJ TTU! 2tTUJ 

sm 1 sin 



(50) 

(51) 



where v*^ = ^{0){1 — cosOy. This rate can be also ex- 
pressed through the angular harmonics Vk of the elastic 
scattering rate v{0) as u*^ = (Si^o — 4z/i -|- 1/2) /2. 

The phonon-assisted (index ph) contribution comprises 
three terms. The first is described as the phonon-assisted 
resistivity modified by heating: 



2m 



-2c;"7^AQ(0)sin 



• (52) 



T . Finally, the third phonon-assisted term can be at- 
tributed to the displacement mechanism of MW photore- 
sistance: 



dis 
Pph 



d P^S2 < J"aq(-w)cos 



2n{iO\Q — cj) 



+>^xq(w) cos 27aq(0) cos — 

U)c U)c 

-^aq(-w) sm ^ H\q{uj) 

. 27r(wAQ + w) . 27rWAQ \ 

X sm — ^ — h 2ftAQ(0) sm — > . (56) 

U!c l^c J 

In contrast to the impurity-assisted contribution pf^, 
this term contains oscillating functions of (waq ioj)/u)c 
describing interference of MIRO and PIRO. 

The next step is to consider the terms associated with 
the oscillating part of the distribution function, Sf^, gen- 
erated by microwaves. Taking into account that Sf^ is 
proportional to MW power, one should retain only the 
terms with n = in Eq. (29). As a result. 



6p. 



m 



deSRl"" + 2Si 



deSRPj' 



,(57) 



where SR^J" and SR^^ are obtained by a linearization of 
R^Jq and R^q with respect to 6/^: 



5RI 



D. 



2d5fe 



de 



(58) 



Here and below, 



1), 



(53) 



^AQ(a;) = (iV..,-iV^,,,J?^%±^. (54) 

The term does not oscillate as a function of w. In 
the limit Te = T it is reduced to the resistance pph cal- 
culated in Ref. 18. The part of p"^ proportional to the 
squared Dingle factor describes PIRO modified by heat- 
ing (note that ?/aq(0) is nonzero only at Te 7^ T). The 
second phonon-assisted term is the MW-induccd "clas- 
sical" photoresistance, which is not related to quantum 
oscillations of the density of states and, therefore, does 
not oscillate with the magnetic field: 



Pph 



2e^n, 



■S2 {J'xQioj) + J-aq(-w) - 2J-aq(0)} . (55) 



This term essentially requires inelastic scattering: a for- 
mal limiting transition waq makes p^f^ equal to zero. 
Therefore, p^f^ has no analogue in the impurity-assisted 

tansport. The contribution (55) is strongly suppressed 
with decreasing frequency w and increasing temperature 



6Rf = DeD, 



dSJe 



x[(iV^ 



-iN^xc+.f:)Sfe^uxc]- (59) 

The integration over energy in Eq. (57) is done with the 
aid of Eq. (44) for 5fs and under the approximation (48). 
The result for dp appears to be rather complicated. It is 
presented below in a simplified form omitting the terms 
quadratic in electron-phonon scattering contributions: 

r "^0^2 tr r„ 27ra; . 27ra; 



^^'^zm— Sin — 



, .C I 27rWAQ . 27rWAQ , „ 

+4So <! — sm —{N^^Q 



UJc 



+2P^,5i ?{AQ(-w)sin 



27r(aJAQ - 1^) 



OJc 

27r(wAQ + w) . 2-KUJxQ 

+^aq(w) sm — ^ — ^ - 2Waq(0) sm ^ 

X cos^ — + (2J-aq(0) - Txdi-uj) - -Faq(w)) 
xcos cos ^ + o (•^aq(w) - ^aq(-w)) 



. 2-KU) . 27rWAQ 1 

X sm sm — > 



.(60) 
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This form is valid under the condition vj^ ^ WpW (see 
the definition of the rate Vpj^ below and the correspond- 
ing discussion in the next section). The first term in 
the square brackets of Eq. (60) is the impurity-assisted 
scattering contribution.^ The second one is the phonon- 
assisted scattering contribution which is caused by heat- 
ing and docs not oscillate with MW frequency. The third 
one (proportional to Pui) is the phonon- assisted scat- 
tering contribution showing interference of MIRO and 
PIRO. Since only the terms linear in MW power are con- 
sidered, the functions and H in Eq. (60) should be 
taken at equilibrium, = T. The same statement con- 
cerns Eq. (56). 

The results are considerably simplified under the con- 
dition 2T 3> w. In this case one can write a general 
expression for the whole resistivity (47), which includes 
also the terms omitted previously in Eq. (60). It is con- 
venient to represent the resistivity as a sum 



p = pb + pos, 
where the background (non-oscillating) part is 



m 

Pb = ^— 



(61) 



(62) 



by the second and the third terms in Eq. (63), which 
are proportional to P^. These terms express contribu- 
tions of the displacement and inelastic mechanisms of 
photoresistance, respectively. Their frequency depen- 
dence is given by the same oscillating functions as in 
the case of impurity-assisted transport, while the 
scattering rates are modified by electron-phonon inter- 
action: v*^ v*^ + v^l and ->■ v^L + v^l. This 



modification is very essential because the rates v^l and 
oscillate with the magnetic field. The interference 



of phonon-induced magnetoresistance oscillations with 
MIRO is produced as a result of multiplication of these 
rates by the frequency-dependent oscillating functions. 
Finally, the last term in Eq. (63) has no analogue in 
the theory of impurity- assisted magnctotransport. It is 
proportional to — T and can be described as a re- 
sult of MW heating on phonon-assisted magnetoresis- 
tance. This term does not oscillate with frequency lo 
but it oscillates with the magnetic field, similar to the 
equilibrium magnetoresistance. The phase of these oscil- 
lations is shifted with respect to the phase of equilibrium 
PIRO. Notice that the heating factor Te/T - 1 is ex- 
pressed through R_j with the aid of Eq. (41), so the whole 
MW-induced correction to p is proportional to P^. 



and the oscillating part includes four terms: 



Pos — 2 



TTOJ TTUJ 



271 UJ 



sm 



— sm 



, tr c1n2 27'"W . 2mj 
-Pu,Tin{yZn + i^ph) — Sin — 



T 



Hl + 2y%nn{yZn + Kh)\\- (63) 



The phonon-assisted transport rate i'*^ is given by Eq. 
(43), while the oscillating rates v^f^ and are defined 



as 



cn _ <j 



r 2T 



■ cos ■ 



AQ 



(■ 



2T 



(64) 



As seen from this definition, oscillations of and 
with the magnetic field have different phases governed 
by cosine (index c) and sine (index s) functions, respec- 
tively. The rate differs from i/*^ by the presence 
of the oscillating factor cos(27ra;;^Q/wc) under the inte- 
gral operator. This rate describes equilibrium PIRO, the 
first term in Eq. (63), in the approximation of over- 
lapping Landau levels {(P <C 1) used here (notice that 
in Ref. 18 v'zl, is denoted as v. 



Xp^j ic uv^iiubtju aa ix^^ ). The dircct influ- 
ence of microwaves on magnetoresistance is described 



IV. NUMERICAL RESULTS AND DISCUSSION 

As follows from the above consideration, quantum os- 
cillations of resistivity in 2D electron layers irradiated by 
microwaves demonstrate a rich and complicated behav- 
ior if electron-phonon interaction is taken into account. 
The absorption of microwaves in the presence of this in- 
teraction leads to the appearance of various additional 
oscillating contributions to magnetoresistance, Eq. (63). 
These contributions are similar to equilibrium PIRO in 
the sense that their periodicity depends on characteristic 
phonon frequencies lo\q with Q = 2pp. However, they 
are essentially non-equilibrium and some of them (those 
entering the last term of Eq. (63)) have a phase differ- 
ent from that of equilibrium PIRO. The non-equilibrium 
phonon-induced contributions entering the second and 
the third term of Eq. (63) show the periodicity deter- 
mined by combined frequencies such as w ± waq and cor- 
respond to the interference of PIRO and MIRO. 

To illustrate the influence of microwaves on magnetore- 
sistance and to compare contributions of different terms 
in Eq. (63), numerical calculations have been carried out. 
The sample parameters are taken for [001]-grown GaAs 
quantum wells of high mobility (1.17 x 10^ cm^/V s at 
low temperature) experimentally investigated in Ref. 17 
without MW excitation. The list of parameters neces- 
sary for description of electron-phonon interaction can 
be found in Ref. 18. The inelastic scattering time was 
estimated according to the relation* 



Tin 



(66) 
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where Ai„ is a numerical constant of order unity. For cal- 
culations, Xin = 1 was chosen, which is consistent with 
experimental data.^^ The rate i'*„^ characterizing the dis- 
placement mechanism of photoresistance under impurity- 
assisted scattering has been calculated with the aid of 
the model of long-range random impurity potential de- 
scribed by the correlator w{q) cx exp(— Zc?)- The cor- 
relation length, Ic ^ Pp^, is determined by using the 
transport rate i/f^ found from the given mobility and 
quantum lifetime — 15 ps estimated in Ref. 17. 
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FIG. 1: (Color online) Upper panel. Magnetic-field depen- 
dence of the oscillating part of resistivity, pos, calculated for 
the quantum well of Ref. 17 at T = 6 K in equilibrium (black, 
curve 1) and under 135 GHz MW excitation with E^j — 2 
V/cm (blue, curve 2) and Eu, = 3 V/cm (red, curve 3). Lower 
panel. MW-induced contributions to pos ett Eu, = 3 V/cm: 
displacement mechanism, second term in Eq. 63 (green, curve 

1) , inelastic mechanism, third term in Eq. 63 (magenta, curve 

2) , non-equilibrium PIRO, fourth term in Eq. 63 (orange, 
curve 3). The contribution due to elastic impurity scattering 
(MIRO) is shown by black dashed line. 

Figure 1 shows the oscillating part of resistivity at a 
lattice temperature of 6 K without microwaves (in equi- 
librium) and under 135 GHz excitation whose intensity 
is determined by the MW electric field E^^. For the 
chosen frequency, the increase in intensity leads to en- 



hancement of PIRO maxima at 0.12 and 0.18 T, while 
the maximum at 0.3 T is suppressed and eventually in- 
verted, and two additional maxima appear at its sides. 
The inversion of this PIRO peak cannot be explained as 
a result of simple superposition of phonon-induced and 
MW-induced oscillations and is attributed to interference 
of these kinds of oscillations, as refiected in Eq. (63). The 
lower panel of Fig. 1 shows relative contributions to the 
oscillating resistivity. Apart from the impurity-induced 
contribution (MIRO), this plot shows the total contri- 
bution of displacement mechanism, inelastic mechanism, 
and non-equilibrium PIRO. All these contributions are 
comparable, though the inelastic mechanism contribu- 
tion is smaller than the others in the region around 0.3 T. 
The i?-dependence of this contribution closely follows the 
MIRO plot, showing slight deviations caused by phonon- 
induced oscillations. In contrast, such deviations are very 
pronounced in the displacement mechanism contribution 
(see, for example, peak at 0.25 T and minimum at 0.4 
T). The contribution of non-equilibrium PIRO is clearly 
shifted by phase with respect to the other contributions. 
This shift leads to a partial compensation of different 
contributions, which reduces the influence of microwaves 
on magnetoresistance. 

The effect of microwaves of a smaller frequency, 50 
GHz, is illustrated in Fig. 2. The magnetoresistance is 
considerably changed only in the low-field region. These 
changes are similar to those in Fig. 1: PIRO peaks 
are either enhanced or inverted by microwaves. Again, 
the inelastic mechanism contribution is weakly modified 
by electron-phonon interaction, while the displacement 
mechanism contribution is strongly modified by this in- 
teraction and shows numerous oscillations related to in- 
terference of PIRO with MIRO. The non-equilibrium 
PIRO contribution is smaller than the others. The rela- 
tive decrease of this contribution with lowering frequency 
also follows from Eqs. (63) and (41). 

The calculations demonstrate that the displacement 
mechanism of MW photoresistance is much stronger in- 
fiuenced by electron-phonon interaction than the inelas- 
tic mechanism. This property is understood from Eq. 
(63) showing that the relative contributions caused by 
phonons for these two mechanisms are determined by the 
ratios and respectively. At T = 6 K 

the second ratio is small, as also seen from experimen- 
tal data:^^ the amplitudes of equilibrium PIRO governed 
by the rate are small compared to the background 
resistivity. On the other hand, for long-range impurity 
potentials one has v*^ ^ t'*^, while v^f^ is approximately 
twice larger by amplitude than i/^j^ , since the main contri- 
bution to the oscillating phonon- assisted rates zy^JJ comes 
from backscattering processes, ~ tt. Therefore, one 
always has 



> 



(67) 



With increasing temperature, when electron-phonon 
scattering gets stronger and the rates z^^JJ; increase, the 
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FIG. 2: (Color online) The same as in Fig. 1 for 50 GHz MW 
excitation with Eui — 0.5 V/cm (blue, curve 2) and Eui — 0.8 
V/cm (red, curve 3). The MW-induced contributions in the 
lower panel are shown for E^i =0.8 V/cm. 



ratio \l'pf^\/l'*^ becomes comparable to miity and the dis- 
placement mechanism contribution is considerably mod- 
ified by this scattering. This corresponds to the tem- 
perature chosen for the calculations. Further increase in 
temperature leads to an interesting situation when the 
amplitudes of v^f^ become much larger than i^*^^, so the 
displacement mechanism of MW photoresistance is dom- 
inated by electron-phonon scattering. 

It is worth noting that at low temperatures (T ~ 1 
K), when phonon-assisted contributions are yet frozen 
out, the inelastic mechanism is parametrically stronger^ 
than the displacement one, because of large value of 
Tin- However, in the samples of very high mobility, 
like those investigated in Ref. 17, the contribution of 
the displacement mechanism experimentally proves to 
be important^^ starting from T = 2 K. According to 
estimates, at this temperature the product Tinvf^ de- 
termining the strength of inelastic mechanism is smaller 
than unity. As follows from Eq. (63), the parameter 
Tinvfm ^Iso determines the relative strength of phonon- 
assisted contribution in the inelastic (third) term com- 
pared to the similar contribution in the displacement 



(second) term. Therefore, in the temperature region 
when phonon-induced oscillating contributions to resis- 
tivity become thermally activated (T > 2 K), their rel- 
ative influence on the MW photoresistance via inelastic 
mechanism is already suppressed. In particular, at T = 6 
K a strong inequality. 



(68) 



allows one to neglect the phonon-induced oscillating pho- 
toresistance due to inelastic mechanism. For the samples 
of moderate mobility (~ 10® cm^/V s), the applicabil- 
ity of Eq. (68) requires higher temperatures. Neverthe- 
less, in the low-temperature region {Tinvlm > 1) the ratio 
Wph\/^im remains small, and phonon-induced photore- 
sistance due to inelastic mechanism can be neglected in 
comparison with the impurity-induced photoresistance. 
Notice also that the condition (68), being applied to non- 
equilibrium PIRO (the last term in Eq. (63)), allows one 
to neglect the part which comes from Sp and is propor- 
tional to Tin- Thus, under condition (68) one may always 
ignore the contribution to Sp caused by electron-phonon 
interaction. In contrast, the non-equilibrium contribu- 
tions to Pe caused by this interaction cannot be ignored 
and are important for evaluation of MW-induced mag- 
netoresistance in 2D systems where PIRO is observed. 
These contributions are described by the oscillating rates 
lyfh and i^'pl in Eq. (63).^ 

The main approximations used in the paper are dis- 
cussed below in more detail. 

The approximation of overlapping Landau levels corre- 
sponds to a simplified (single-mode) description of the os- 
cillating density of electron states in contrast to the exact 
SCBA expression (25). Though in high-mobility samples 
separation of Landau levels becomes essential already at 
i? ~ 0.1 T, this description produces a very good agree- 
ment with the calculation based on Eq. (25) (see Ref. 
18) for equilibrium oscillating magnetoresistance up to 
B = 0.4 T. This means that the first PIRO harmonic 
dominates in the main interval of magnetic fields under 
consideration, and only the last PIRO peak at _B ~ 0.6 
T appears to be slightly sensitive to the shape of the 
density of states. Besides, regarding a limited validity of 
the SCBA in the regime of separated Landau levels, the 
single-mode description proves to be a reasonable choice. 
It also provides a rigorous justification for introduction 
of inelastic scattering time necessary for analytical de- 
scription of inelastic mechanism of photoresistance. 

The assumption of weak MW excitation means that 
only linear terms in the expansion of magnetoresistance 
in powers of MW intensity are taken into account. First, 
this implies a neglect of multi-photon absorption pro- 
cesses and formally requires that the argument (3 of the 
Bessel function in Eq. (29) must be small for the scatter- 
ing angles essential in the integration. Since the phonon- 
assisted scattering occurs in a broad range of angles, the 
condition /3 ^ 1 is equivalent to P^j <C 1 and is satisfied 
by a proper choice of the MW field strength E^- Next, 
the neglect of saturation effect for inelastic mechanism 
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contribution*'^ requires Tinvl'^P^ <Si 1. As shown above, 
the phonon-induced photorcsistance oscihations in high- 
mobiUty layers are essential at Tinvl^^ < 1, so the neglect 
of saturation is also satisfied at P^^ ^ 1. Finally, the 
absence of strong heating, Te/T—1 <^ 1, is equivalent to 
smallness of the right-hand side of Eq. (41) and also can 
be expressed in terms of smallness of E^^ . 

Coming back to multi-photon processes, it is worth 
noting that phonon-assistcd multi-photon absorption re- 
quires much weaker MW fields than impurity-assisted 
multi-photon absorption. The reason is the long- 
range nature of impurity potential in high-mobility 
(modulation-doped) structures, which leads to small- 
angle scattering (^ <C 1) and, thereby, effectively smaller 
/3 in the impurity term of Eq. (29). Since the multi- 
photon processes are thought to be responsible for the 
phenomenon of fractional MIRO observed at elevated 
MW intensities, the above conclusion may be 
useful for interpretation of relevant experimental data. 
In particular, a recent experiment'''* shows several high- 
order fractional resonances in photorcsistance when the 
MW power is still insufficient to produce these features 
as a result of multi-photon processes under electron- 
impurity scattering. A consideration of electron-phonon 
interaction, which enables multi-photon processes at 
weaker MW power, could possibly explain the data of 
Ref. 34 (indeed, the high-order fractional features are 
observed at T = 6.5 K, when phonon-assistcd scattering 
may become important). A more detailed consideration 
of this problem requires a separate study which is beyond 
the scope of the present paper. 

The neglect of electron-phonon and electron-electron 
interaction in calculation of the Green's function in Eq. 
(14) is justified at low temperatures, when the inverse 
quantum lifetime due to impurity scattering, l/xg, is 
much larger than the corresponding inelastic scattering 
rates. This leads to temperature-independent density of 
states characterized by the Dingle factors used in the 
above calculations. However, while electron-phonon con- 
tribiition is not essential in a wide tcmperatiire range, the 
contribution of electron-electron scattering to the Landau 
level broadening cannot be disregarded already at several 
Kelvin. To take this effect into account, one should re- 
place the Dingle factors d with temperature-dependent 
functions d{T) = exp[—'jT/ucTq{T)], where 



1 



AT) 



1 1 
— + — 



(69) 



and the numerical constant A (usually determined exper- 
imentally) is of the order of unity. The reliability of this 
approach is justified theoretically and proved in numer- 
ous works by studying temperature dependence of quan- 
tum oscillations in magnctorcsistancc.*^''^^'^*'^^'^^'^*'^^ 
The inclusion of temperature-dependent Dingle factors 
into Eq. (63) leads, in particular, to non-monotonic tem- 
perature dependence of PIRO amplitudes, which is also 
observed in experiments. ^^'^"^ 



The consideration was done for the case of spatially 
homogeneous 2D electron system. The effects of a finite 
size of the sample, in particular, the infiucncc of con- 
tacts and edges on the MW photorcsistance, have been 
neglected. This includes, for instance, the neglect of 
confined magnetoplasmon effects in MW absorption. ^-"^ 
The problem of magnetoresistance under inhomogeneous 
MW field distribution appears to be even more impor- 
tant. The observed inscnsitivity of MIRO to the sense of 
circular polarization,'*^ the absence of MIRO in contact- 
less measurenents,^^ and calculations of the field distribu- 
tion in 2D systems with metallic contacts^'' suggest pos- 
sible near-contact origin of MIRO phenomenon in high- 
mobility layers. This means that, owing to large MW 
field gradients in the vicinity of contacts, there may exist 
an additional (caused by near-contact regions) oscillating 
contribution to resistivity*^ which exceeds the contribu- 
tion caused by uniform MW field in the bulk of thc^ layer. 
The nature of the assumed near-contact oscillating con- 
tribution remains unclear, and its theoretical description 
is currently missing. Further investigation of this prob- 
lem is necessary. In the meanwhile, any extension of the- 
oretical knowledge about MW-induced response based on 
bulk properties of 2D electron gas is important and may 
help to uncover the exact nature of MIRO. 

In summary, a theoretical study of the effects of con- 
tinuous MW irradiation on the magnetoresistance of 
2D electrons interacting with impurities and acoustic 
phonons is presented. The main subjcict of the study 
was the quantum oscillations of resistivity caused by the 
scattering-assisted electron transitions between different 
Landau levels. In the absence of microwaves (in equilib- 
rium) , there is only one kind of such oscillations (PIRO) 
owing to electron-phonon scattering. When MW radia- 
tion is applied, the transitions between Landau levels can 
occur with absorption or emission of MW quanta, either 
under elastic scattering by impurities or under inelastic 
scattering by acoustic phonons. The concept of micro- 
scopic mechanisms of MW photorcsistance (displacement 
and inelastic mechanisms), previously developed for the 
case of electrons interacting with impurities, applies to 
both these kinds of transitions. The impm-ity- assisted 
transitions arc responsible for the oscillations (MIRO) 
whose 1/B periodicity is determined solely by the MW 
frequency u. The phonon- assisted transitions lead to 
more complicated oscillations which arc described as a re- 
sult of PIRO-MIRO interference. It is demonstrated that 
such oscillations are attributed mostly to the displace- 
ment mechanism. Finally, heating of electron system by 
microwaves gives rise to an additional phonon-induced 
oscillating contribution, denoted here as non-equilibrium 
PIRO, which is shifted by phase with respect to equilib- 
rium PIRO. This contribution is not related to specific 
features of MW excitation and appears also when elec- 
tron system is heated by a dc field. Accounting for all 
the oscillating contributions leads to a peculiar magne- 
toresistance picture that should be observable in high- 
mobility 2D layers if the temperature is high enough so 
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the electron-phonon scattering is not frozen out. Some 
relevant examples are given in Figs. 1 and 2. 

One of the primary goals of this work was to em- 
phasize the importance of electron-phonon interaction 
in description of the microwave-induced oscillating 
magnetoresistance. The author believes that the results 
and conclusions of this research will be useful for better 



understanding of existing experimental data and may 
stimulate further experiments and theoretical studies. 

Acknowledgement: The author is grateful to I. A. 
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